Quantum mechanical inverse scattering problem 
at fixed energy: a constructive method 



t 



Tamas Palmai* and Barnabas Apagyi 



Department of Theoretical Physics 
Budapest University of Technology and Economics 
Budafoki ut 8., H-1111 Budapest, Hungary 

November 28, 2011 



Abstract 

The inverse scattering problem of the three-dimensional Schrodinger 
equation is considered at fixed scattering energy with spherically symmet- 
ric potentials. The phase shifts determine the potential therefore a con- 
structive scheme for recovering the scattering potential from a finite set of 
phase shifts at a fixed energy is of interest. Such a scheme is suggested by 
Cox and Thompson [3] and their method is revisited here. Also some new 
results are added arising from investigation of asymptotics of potentials 
and concerning statistics of colliding particles. A condition is given [2] for 
the construction of potentials belonging to the class Li.i which are the 
physically meaningful ones. An uniqueness theorem is obtained [2] in the 
special case of one given phase shift by applying the previous condition. 
It is shown that if only one phase shift is specified for the inversion proce- 
dure the unique potential obtained by the Cox- Thompson scheme yields 
the one specified phase shift while the others are small in a certain sense. 
The case of two given phase shifts is also discussed by numerical treat- 
ment and synthetic examples are given to illustrate the results. Besides 
the new results this contribution provides a systematic treatment of the 
CT method. 
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1 General results, introduction 

We start with the Schrodinger equation in at a fixed positive energy 

in R3 (1) 



x\ oo, a' = — , (2) 
r 

where a is the direction of the incident wave and A{a,a') is the scattering 
amphtude. The following general theorem is due to Ramm |8j. 

Theorem 1.1. A(a',a) scattering amplitude \fa' € Sf, Va € S'| (arbitrary 
small open subsets of S'^) determine q{x) uniquely in the function class Qa = 
{q:q = q, q{x), \x\ > a, q{x) £ L'^{Ba)}, Ba = {x : x e M^ < a}. 

The present treatment is restricted to spherically symmetric potentials, that 
is q{x) = q{r). In this case we have the partial wave expansion of the wave 
function: 

oo £ j { \ 

*(x,a)=^ 47ri^^^li™(x/r)y,„(a), (3) 

where il}i{x) satisfies the radial Schrodinger equation (see below) with the ap- 
propriate boundary conditions. 

For spherically symmetric potentials the scattering amplitude takes a similar 
expansion form, namely 

oo t 

A{a',a)= A{a' ■a) = Y^ ^ AtYt^{a')Yt„,{a), 

The phase shifts, {^£}f=o,i,2,... can be restricted to — ^ < < f , since from 
the experimental point of view they are undetermined to an additive factor 
of kn, fc G Z. Then it is apparent that the knowledge of the phase shifts is 
equivalent to that of the scattering amplitude. However not all the phase shifts 
(as considered in the theorem of Loeffel pLi) are needed for unique reconstruction 
of the potential. The following result is due to Horvath and Ramm [9l |6] . 

Theorem 1.2. The phase shifts {(S^j^go determine the potential uniquely in the 
class Q — {q ■ q ~ q, q{x) — q{r), r = \x\, q{r) — 0, r > a; rq{r) G Li(0,a)} if 
the Miintz condition holds for £: 

E ! = (4) 

£££.,19^0 

This condition is almost necessary in the sense that in the class Qa with < cr < 
2 the set {5t} with ^^"^ < oo is not enough to recover the potential uniquely. 
Qa = {q ■■ q ^ q, q{x) = q{r), r ^ \x\, q{r) = 0, r > a; r^^'^q{r) G Li(0, a)}. 



[V^ + 1 - (7(a;)]*(.T,a) = 
whose scattering solution takes the form 

e*"" / 1 ' 

^(x,a) = e"" + A(a',a) — +o - 

r \ r 
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It is because of this result that one may look for a potential when a part of 
the scattering phase shifts are known. 

The phase shift 6e appears in the asymptotic form of the regular solution of 
the Schrodinger equation, which is defined in the following manner: 



T + T 



(2^+1)!! 



0, 



in 



ipi{r) = Bi sin(r - y + ^f) + o(l), 



r — > oo 



with B( being a constant. Note that for the q = zero potential we have 

LroUi{r) = i{£ + l)ui{r), where L^o = r'^-^+r'^ and Ui{r) = ^^J^+i (r) 

and J„(r) is the nth order Bessel function of the first kind. 

The Povzner-Levitan representation in the fixed energy problem [5] has long 
been used, however the existence and uniqueness of a K{r,r') £ C^(R+ x 
transformation kernel have been proven rigorously only recently for potential 
q{r) G C^(0,a) independent of £ and satisfying 



ipe{r) = ui{r) 



(5) 



For K{r, r') we have a Goursat-type problem (equivalent to the Schrodinger 
equation) : 

LrK{r, r') = Lr'oK{r, r'), < r' < r, 
2 d K{r,r) 



q{r) 



r dr 



K{r,0) = 0. 



Now, in analogy with the fixed-£ problem a Gel'fand-Levitan-Marchenko- 
type integral equation is written up for K{r, r'): 



K{r,r') = g{r,r') - / dpp '^K{r, p)g{p,r'), r > r' . 



(6) 



Since K{r, r') is unique this integral equation must have a unique solution 
for K{r, r') to make sense. Furthermore, for consistency with the original 
Schrodinger equation, g{r, r') must satisfy 



L,off(r, r') - Lr'^gir, r'), g{r, 0) = 5(0, r') = 
and through the integral equation 

2d K{r,r) 

w) = — 

r dr r 



(7) 



(8) 
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is maintained. g(r,r') is a good candidate for approximation since its differen- 
tial equation can be satisfied trivially. For instance if the angular momentum 
expansion ^5, 

9{r,r') ^^cae{r,r'), (9) 

e 

is imposed then for the 'ye{r, r') functions we have only the potential-independent 
restrictions of 

Lr^li{ry)^Lr.oli{ry), 7,(r, 0) = 7,(0, r') = 0, VI (10) 
and the information is only stored in the ci expansion coefficients. 

2 Cox— Thompson (CT) method 

In the framework of the method proposed by Cox and Thompson ^3, one takes 
the following separable form for the 7,(r, r') functions 

r') = u^(min(r, r'))u£(max(r, r')), (11) 

which is also the Green's function of the q = {) radial Schrodinger equation for 
the ^th partial wave: 

' <f . m + iv 



li{r, r') = 6{r ~ r'); (12) 



and the summation in g{r,r') runs only over a finite set S of ^'s: 

gi.'T^ '"') = X! C£M£(min(r, r'))ve{max{r, r')), (13) 

containing the Riccati-Bessel functions, connected to the Bessel and Neumann 
functions by 

Mr) - \I^JiH^r), v,{r) ^ ^J^Y,^.{r). (14) 
For solving the GLM-type equation ^ we use the separable ansatz 

Kir,r')^J2^L{r)uLir') (15) 
LeT 

with a finite set T of "shifted angular momenta" satisfying S' n T = and 
|S'| = \T\. The use of such an ansatz can be motivated by the following result 
ascertained from [4j. 

Proposition 2.1. If r'^^^^K{r,r') = 0(1), r' ^ is imposed on K{r,r') and 
the L numbers are restricted to L > —0.5 then 

llL<,TW+l)-LiL + l)) 



^ K{r,r') = AL{r)uL{r'). (16) 



m'es,£'#£W^+l)-^'(^' + l)) 
Where {c,} O {L} is a one-to-one mapping. 
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Also, this form provides a reasonably easy way to solve the GLM-type inte- 
gral equation. In fact it makes a system of algebraic equations instead of the 
integral equation, namely 

The parameters of the set T can be obtained from the phase shifts {Si}g^s 
through the transformation equation ([5]) which in terms of the L's takes the 
form 

■^-^ , , ^UL(r)u'i,(r) — u't (r)ui(r) , ^ 

Mr) = - E^.W l^i:i)-LiLil) 

Taking both equation fUEh and pT)) for large r's, i.e. at r — s> oo we get the 
system of nonlinear equations |TT1[T^[T3] connecting {Sg}t^s to T: 

^2..,^!+^^ (19) 
1 — i/C„ 



with 



= E E [Msin]iL[M-,l]u'e^^'^'-'>^' e G S, (20) 
Msin \ _ 1 / sin((£-L)7r/2) , K^SLeT. 



M,os J L{L + 1) - ^(^ + 1) \ cos {{£ ~ L)ti/2) 

(21) 

After some investigation of this nonlinear system of equations one can conclude 
that for a given set of phase shifts it yields an infinity of solutions, each giving rise 
to a potential. It was our result in 2^ that one can select a physical potential 
(perhaps uniquely) with the aid of a consistency check. This check is first 
revisited in Sec. 4 and then it is applied for the special cases of one and two 
dimensions, \S\ = \T\ = 1, 2. 



3 Some new results concerning the CT potential 

Before continuing with the consistency check it is worthwhile to further study 
the CT method particularly concerning the potential one can obtain by applying 
it. We address the case when only a finite number of input phase shifts are used 
to construct the potential. 



3.1 Asymptotics 

First, we show that the potential is generally not compactly supported nor is of 
long-range. To see this define the functions {ylJ^(r)}LeT by the limit 

Adr)^ Al{r)^-o{l), r ^ oo, LeT. (22) 
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We infer from that the asymptotic coefficients A1 (r) assume the form 

AUr) = ql COST + bLs'mr (23) 



where a^'s and 6l's are constant depending on all the elements of both S and 



cos(f(£-L)) 



-^sinf^^ 
1 V 2 



e S. 



(24) 
(25) 



This form in turn implies that 

i4:(r, r) = asin(2r) +/3cos(2r) +7 + o(l), r oo (26) 

and thus 



, , , /3sin(2r) — acos(2r) , 
q{r) = ^ + o(r-2), 



r — ?> oo 



(27) 



which means that the potential generally falls of like an inverse power of two. 

One can give a necessary condition for the potential to decrease more rapidly 
then 0{r~^). One only needs 

a = ^ = 0. (28) 
The quantites a and /? are given by 



a = i ^ {cLL cosL^ - 5l sinL^ 



LeT 



^ = ^ (oLsinL^ + 5lcosL^ 



(29) 
(30) 



Alternatively, K{r, r) can be written as a series in £ (this formula does not hold 
for K{r, r') with r ^ r'), 



K{r,r) = y^^ci(pi{r)vi{r), 



(31) 



ies 



which can readily be seen from the GLM equation (|6]) taken at r = r' and the 
CT formula (fT3|) substituted for g{r, r'): 

K{r,r) ^^^CiUi{r)vi{r) - dpp^^ K{r, p)'^ciui{p)v£{r) (32) 



ui{r) - / K{r,p)ui{p)p "^dp 



(33) 



where the formula ([5]) for (feir) has appeared. This allows for the alternative 
conditions 

'^{-lYaBeCOsSi^O, and ^(-1)^^^^ sin^^ = (34) 
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involving the expansion coefficients, the input phase shifts and the normalization 
constants of the partial wave functions. 

These new results may serve as useful tools to check numerical results or 
incorporated into a solution method they might provide a way to control the 
undesirable oscillations of the inverse potential. 

At this point we shortly discuss the applicability of the CT method. The 
inverse potential is finite at the origin starting with a zero derivative (see e.g. 

and possesses also a finite first moment rq{r)dr < oo. At fixed scat- 
tering energy the CT procedure is thus a particularly successful method for 
(re)constructing inverse potentials. It was used to recover interaction potentials 
(see e.g. [TT1[TU[T3]) for various physical system. Furthermore it is also possi- 
ble to generalize the method to recover potentials having Coulomb tail e.g. by 
replacing the Riccati-Bessel function by Coulomb wave functions [15l [16] . 



3.2 CT potentials regarding statistics of colliding particles 

In real-life scattering experiments one encounters cases when the colliding par- 
ticles are zero spin bosons. Then it is well-known that only the even numbered 
partial waves contribute to the scattering amplitude because of the symmetry 
of the wave function. In such cases with the CT method it is possible to use 
only the experimentally available phase shift data corresponding to even partial 
waves. It is then reasonable to ask for the value of phase shifts provided by the 
CT potential for the odd partial waves. It turns out that these odd phase shifts 
are zero. 

To see this let us derive the phase shifts of the CT potential which can be 
obtatined from Eq. (fT8)) taken at r — >■ oo: 



Now the trigonometric form ((23)) is substituted for (r) . There are two special 
cases that we consider: 

i) the set S of input angular momenta consists of only even numbers, 

ii) S consists of odd numbers. 

In the first case we have — Q VL G T while in the second = ML e T. 
This is deductible from Eqs. p4)) and ([25]) which assume the forms 



E 

LeT 

E 



aL _ 

bL j L{L + 1) - £{£ + 1) ^ \ 

ai \ sin (if) _ f 



bL j L{L + 1) - £{£ + 1) \ 1 



£ e S (36) 
£eS (37) 



for the cases i) and ii) respectively. One can assume cos (if) 7^ and sin (L-^) 7^ 
0, since the matrix AI with elements M^l — {L{L + 1) — £{£ + 1))^^ is invertible 
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(it is a Cauchy matrix) and cannot be singular unless T C\ S ^ % which is ex- 
cluded by assumption. We thus have 6l = VL £ T for i) and = ML £ T 
for ii). Now 



which implies 

fj{r) — Bj sin(r— j7r/2) + o(l) for i) and odd j, or for ii) and even j. (39) 

In other words the CT phase shifts of the opposite parities are exactly zero. 
Notice that the CT method allows the construction of potentials which are 
transparent for half the partial waves (being even or odd in parity) . 

Applications suggest that if dealing with input partial wave data of one parity 
the performance of the CT method with the same number of input diminishes 
compared to the case when data with both parity is employed. Therefore our 
sum rules derived in the previous subsection can be extremely useful to improve 
performance by supressing oscillations of the potential. Also note that one 
of our sum rules simplify (due to Bi cos 5^ — I, £ £ 5 in the even case and 
Bi sin Si = 1, ^ G S* in the odd case [13]) to 



while the other becomes '^^^g cgtandi = and '^^^g ci cot Si — for the even 
and the odd case, respectively. 

4 Consistency check 
4.1 General condition 

The kernel g{r, r') only makes sense if the integral equation is uniquely solvable 
with it. As Eq. ([6]) can be viewed as a Fredholm type integral equation of the 
second kind for fixed r, viz. if rr'K{r,r') = K{r,r') and rr'j{r,r') — g{r,r') we 
have 



According to Fredholm's alternative the Fredholm determinant thereof must be 
nonzero for all fixed r > 0. When using the CT ansatz the Fredholm determinant 
of the integral equation becomes the determinant of the system of the algebraic 
equations ([TT)) . 




(38) 




(40) 




(41) 




(42) 



It is easy to see that foi X eft = {x : D{x) = 0, x e K+} one gets 




(43) 
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thus the potential is not in — {q : t\q{t)\dt < oo}. While for = we 
have [3] 

One can conclude the following [5]. 

Theorem 4.1. // and on^j/ if D{r) ^ on r > we get a unique solution of the 
GLM-type integral equation in C^(M^ x M+) and from that an inverse potential 
in Li i = {q : /p t\q{t)\dt < oo}. 

From jjj we know that D{r) 7^ on r > is not the case for arbitrary choice 
of 5* and T: e.g. there it was shown, that if 5 = {0} and L — {2} then D{r) = 
at some r > 0. 

With the help of Theorem 4.1 we can convince ourselves that a particular 
CT inverse potential obtained numerically from arbitrary data is integrable or 
not. If it is so then that potential will generate the input phase shifts. 

Also, using the theorem it is possible to determine the admissible set of L 
numbers for given ^'s. In the one-^ case we get a straightforward admissible 
set of i's, however in higher dimensions the formulation becomes extremely 
involved. 



4.2 One dimensional case 

In this case - S — {£}, \S\ — 1 - the function W{uL,V£){r) = UL{r)v'g{r) — 
u'i^{r)vi{r) must be examined carefully. This was performed in [2] and thus it 
shall not be discussed here in detail. We only give the key idea of the proof: 
the Wronskian 

ULir)v'i,{r) -u'^{r)vi{r) 

e{i + 1) - L{L + 1) ^ ' 

is nonzero on if and only if the constituent Bessel functions { Jl+i/2{x) and 
^^+1/2(2;)) are interlaced. This is deductible from the observation 

e(i+l)-L{L+l) =1 ^^^P^^^P^P ^P- (46) 
We could prove the next theorem of broader interest [TO] 

Theorem 4.2. The positive zeros of the Bessel functions Jy{x), J'y{x), Y,j(x), 
Yl{x), Jiy+e(x), Y^^!;[x) for nonnegative orders, v >Q are interlaced according 
to the inequalities 

V < < 2/1^4 < Vi^+e.l < y'u,i < ju.l < ji^+eS < jl,s+i < • ■ • (47) 

if and only i/O < e < 1 (otherwise Hu+e.s > ju,s for some s and yu,s'+i > ju+e,s' 
for some s' ). 

The following gives the admissible set \t — L\ < 1. 
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Theorem 4.3. W{uL,ve){r) has no roots on r E M^, that is at N = 1 the 
GLM-type equation is uniquely solvable for the CT method with S ~ {£} and 
T — {L} if and only if \L ~ t\ <\. I E (—0.5, oo), L G (—0.5, oo) is supposed. 

This result allows us to choose uniquely from the solutions of the system of 
equations ([T^ as the solution for one input phase shift is 

2 

L^l - -5i + 2k, fc e Z, (48) 

TT 

which for (5^ G [— §, §] in conjunction with the previous theorem yields 

L^£--5i. (49) 

TT 

At low energies it may happen that only one partial wave contributes mostly 
to the scattering amplitude (e.g. for some partial wave resonances). It is worth- 
while to look at the phase shifts yielded by the CT inverse potential at the 
one-phase-shift level to get a sense of the quality of the inversion procedure. 

Deductible from Eq. ([T5|) is the following formula 

fo, ^odd 
iimSi^\ £(£+1) . . (50) 

l£(£+iF?(£+T) tan<5o, £even 

which specifies the phases of the CT potential for = {0}. 

If e.g. the phase shift is restricted to describe an attractive potential (i. e. 
^0 < 0) the bound 

fo, £odd 
tan<5,<<^ 4 , . , (51) 
I tan Oq , (- even 

can be found using Eq. . This result assures the proper reproduction of the 
phase shifts by the CT potential. 

To illustrate these results Fig. 1 shows synthetic test potentials correspond- 
ing to ^ = 0, (5o = 0.27r obtained by the CT method. In addition to the Li i 
potential an inconsistent one is also shown where k in Eq. (j48| is chosen to be 
other than zero. As indicated before we get a non-integrable potential. 

4.3 Two dimensional case 

Already in this case the Fredholm determinant becomes complicated, which is 
also apparent from it's integral representation 

D{r) = / ULAp)uL2{p')ividp)vi2ip') -'Vijp')vi^{p))p-^p'-'^dpdp'. (52) 
Jo Jo 

Therefore instead of the analytical treatment we determine the admissible set 
of i's numerically for some choices of ts. 
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Figure 1: Potentials obtained by the CT method corresponding to £ = 0, 60 = 
0.27r with L = —0.4 (soUd hne) and L = 1.6 (dashed line). 
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■ Integrable 

■ Non-integrable 



1 2 3 4 5 6 7 
r 
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Our numerical method was to check the determinant at the points of a fine 

lattice on the L1-L2 quarter plane of IR+ x IR+ whether it has any zeros on 
(0, A), where A is a large number chosen to be great enough for 



D{r > A) = const. + e(r), 



\e{r)\ < e 



(53) 



with a very small e. This can be done since every D{r) in any dimensions 
IS'I < 00 has only a finite number of zeros, since the constituent Wronskians all 
tend to constants at large r distances, viz. 



W{r) = UL{r)v[{r) — u'i^{r)ve{r) = cos {i — L) 



+ 



r ^ 00. (54) 



On Fig. 2 the admissible sets of {Li, L2} pairs are depicted for the particular 
choices 5i = {1,3} and 5*2 = {1,2}. 

In the next example (Fig. 3) we calculated some possible {Li,L2} pairs for 
a given {Sg-^ , (5^2 }• We note that only one of them is inside the permitted domain 
and could only find a single solution of the system of non-linear equation that 
is permitted by the consistency condition. Again the Li^i and an inconsistent 
potential is shown. 



5 Conclusion 

After recalling some crucial results of the inverse scattering problem for the 
Schrodinger equation at fixed scattering energy we reviewed a particular con- 
structive scheme, the Cox-Thompson (CT) method for recovering scattering 
potentials responsible for a finite set of given phase shifts. 
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Figure 2: The admissible sets of the T elements for (a) Si = {1,3} and (b) 
5*2 = {1,2} denoted by blank areas. 



(a) 



(b) 




Figure 3: Potentials obtained by the CT method corresponding to 5 = 
{0,1}, with phase shifts calculated from a Woods-Saxon potential. T = 
{-0.3056,0.9295} (solid line) and T = {1.0650,1.7016} (dashed line). Also, 
the original potential, q{r) = ~ [l + e^ '^'^''^^^] , yielding the phase shifts 
{So = 0.4389, <5i 0.1246) is depicted (dotted line). 
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We have analyzed the asymptotic properties of the CT potentials and from 
that we have obtained a sum rules for the expansion coefficients (Eq. (|28I34I) '). 
Then we considered CT potentials corresponding to specific quantum statistical 
characters of the colliding pattern. If the scattered particles are zero spin bosons 
(e.g. oxygen atoms or carbon nuclei) the CT potentials are constructed from 
even phase shifts being accessible for measurements. We have shown that in this 
case the CT potentials in addition to reproducing the input even phase shifts 
give exactly zero phases for the odd partial waves. 

Then we turned to the consistency check. It has been shown that for the 
potential to be integrable (or more precisely, to belong to the function class Li.i) 
a condition must be fulfilled for the intermediate quantities L of shifted angular 
momenta of the CT method. We discussed the case when only one input phase 
shift is used in which circumstance we have explicit uniqueness. The two-phase- 
shift case has also been discussed. The condition for the L numbers to produce 
an integrable potentials has been obtained through a numerical calculation for 
two particular cases and the results have been illustrated in Fig. 3. 

References 

[1] J J Loeffel, Ann. Inst. Henri Poincare 8 339-447 (1968) 

[2] T Palmai and B Apagyi, J. Math. Phys. 51 022114 (2010) 

[3] J R Cox and K W Thompson, J. Math. Phys. 11 805-815 (1970) 

[4] J R Cox and K W Thompson, J. Math. Phys. 11 815-817 (1970) 

[5] K Chadan, F C Sabatier, Inverse Problems in Quantum Scattering Theory 
(Springer, 1977) 

[6] A G Ramm, Commun. Math. Phys. 207 231-247 (1999) 

[7] A G Ramm, Applic. Anal. 81 833-836 (2002) 

[8] A G Ramm, Inv. Probl. 4 877-886 (1988) 

[9] M Horvath, Trans. Am. Math. Soc. 358 5161-5177 (2006) 

[10] T Palmai and B Apagyi, J. Math. Anal. Appl. 375 320-322 (2011) 

[11] B Apagyi, Z Harman and W Scheid, J. Phys. A: Math. Gen. 36 4815 (2003) 

[12] O Melchert, W Scheid and B Apagyi, J. Phys. G: Nucl. Part. Phys. 32 849 
(2006) 

[13] T Palmai, M Horvath and B Apagyi, J. Phys. A: Math. Theor. 41 235305 
(2008) 

[14] T Palmai, M Horvath and B Apagyi, Mod. Phys. Lett. B 22 2191-2199 
(2008) 



13 



[15] T Palmai, B Apagyi and W Scheid, J. Phys. G: Nucl. Part. Phys. 37 025101 
(2010) 

[16] D Hilk, W Scheid, T Palmai and B Apagyi, J. Phys. A: Math. Theor. 43 
225302 (2010) 



14 



